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Abstract —In this paper we consider optimal multiuser down¬ 
link beamforming in the presence of a massive number of arbi¬ 
trary quadratic shaping constraints. We combine beamforming 
with full-rate high dimensional real-valued orthogonal space time 
block coding (OSTBC) to increase the number of beamform¬ 
ing weight vectors and associated degrees of freedom in the 
beamformer design. The original multi-constraint beamforming 
problem is converted into a convex optimization problem using 
semideflnite relaxation (SDR) which can be solved efficiently. 
In contrast to conventional (rank-one) beamforming approaches 
in which an optimal beamforming solution can be obtained 
only when the SDR solution (after rank reduction) exhibits the 
rank-one property, in our approach optimality is guaranteed 
when a rank of eight is not exceeded. We show that our 
approach can incorporate up to 79 additional shaping constraints 
for which an optimal beamforming solution is guaranteed as 
compared to a maximum of two additional constraints that 
bound the conventional rank-one downlink beamforming designs. 
Simulation results demonstrate the flexibility of our proposed 
beamformer design. 

Index Terms —Downlink beamforming, shaping constraints, 
semideflnite relaxation (SDR), orthogonal space time block cod¬ 
ing (OSTBC). 

I. Introduction 

With the massive growth of the number of wireless com¬ 
munication users and the increasing demands for high-rate 
services, the spectral resource is becoming more and more 
scarce. Research on spectrally efficient transmission schemes 
for current and next generation cellular networks that are 
capable of mitigating effects of multiuser and co-channel inter¬ 
ference is attracting considerable interest m. As a spectrally 
efficient multi-antenna technique ||2|, downlink beamforming 
has been extensively studied in the past few years 13- 
0. With the aid of channel state information (CSI) at the 
transmitter, downlink beamforming is employed at the base 
station of cellular networks to serve multiple co-channel users 
simultaneously using spatially selective transmission. 

As a pioneering work in downlink beamforming, the authors 
in 13 consider the problem of minimizing the total trans¬ 
mitted power subject to quality of service (QoS) constraints 
in the form of minimum signal to interference plus noise 
ratio (SINR) requirements at each user. A particular form of 
uplink-downlink duality theory is established in 13 and under 
this framework the downlink beamforming problem is solved 
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using a computationally efficient power iteration algorithm. 
A similar approach that exploits uplink-downlink duality is 
proposed in 0, where the downlink beamforming problem 
of maximizing the minimum SINR subject to a total power 
constraint is considered. 

A different class of approaches is presented in a, Q 
and Q where the downlink beamforming problem is ad¬ 
dressed using conic optimization. The authors in a, a 
solve the beamforming problem by resorting to the concept of 
semideflnite relaxation (SDR) and prove that from the rank- 
relaxed problem a rank-one solutiord can always be obtained 
when the problem is feasible. Besides, the authors in 0 cast 
the problem into a computationally efficient standard second 
order cone program (SOCP) and for which the corresponding 
optimality conditions are derived. 

The multiuser downlink beamforming approaches refer¬ 
enced above all optimize the beamforming weights considering 
the SINR requirements of the individual users served in the 
network. In addition to this, supplementary shaping constraints 
on the beamforming weight vectors are embedded in the 
downlink beamforming problem to incorporate a variety of 
requirements in diverse applications il-ES. For example, 
in hierarchical cellular networks operating under the licensed 
shared access (LSA) paradigm ll23 . pico- and femtocell net¬ 
works co-exist in the same frequency band with the surround¬ 
ing macrocell ns-iEi. Shaping constraints are employed at 
the femtocell base stations to limit the power leakage to the 
macrocell users a-ini and the power leakage to concurrent 
femtocell networks US-Ill. Similarly, in the newly emerg¬ 
ing context of physical layer secrecy, shaping constraints are 
applied to guarantee that the SINRs at the eavesdroppers reside 
below a given detection threshold such that the confidential 
information can only be decoded at the desired receiver M- 
El. Recently downlink beamforming has been associated 
with wireless charging in energy harvesting communication 
networks. In this context shaping constraints are employed 
to ensure that the received power at the harvesting nodes 
is greater than a prescribed threshold to facilitate efficient 
charging HS-lEOl. Furthermore, shaping constraints are used 
in multiuser downlink networks to limit the interference power 
leakage to co-channel users, e.g., in neighboring cells ED, 

E3. 

The above mentioned SDR approach lends itself for appli- 

*By a rank-one solution we mean that the solution matrices of the SDR 
problem exhibits the rank-one property. 
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cation in the multi-constraint downlink beamforming problems 
with a large number of additional shaping constraints 141. ||5l. 
However, if the number of additional shaping constraints is 
large, the relaxation is not guaranteed to be tight and a general 
rank SDR solution may be obtained resulting in suboptimal 
or even infeasible beamforming solutions. It is demonstrated 
in m, 123 that when the number of additional shaping 
constraints is upper bounded by two, a rank-one solution 
can always be found by applying a particular rank reduction 
algorithm, without losing the optimality of the solution. In 
other cases, if the solution after rank reduction still exceeds the 
rank of one, a suboptimal solution can be generated from the 
SDR solution by using, e.g., randomization techniques l2^ . 

123. 

In this work, we develop a new approach to optimally 
solve the downlink beamforming problem in the case that the 
number of additional shaping constraints is no greater than 
79. We exploit CSI knowledge at the transmitter and combine 
downlink beamforming with full-rate high dimensional real¬ 
valued OSTBC to increase the degrees of freedom in the 
beamformer design. Several works have proposed the idea of 
combining beamforming with space time coding Ell-ES. 

In l28l . side information in the form of channel estimates 
is utilized to design linear beamformers for OSTBC precoded 
transmission based on a pairwise error probability (PEP) cri¬ 
terion. Two-directional Eigen-beamforming based on channel 
mean feedback is investigated in using beamforming 
along with Alamouti coding 1^ . and the symbol error prob¬ 
ability (SER) criterion is employed in the beamformer design. 
A similar idea is applied in where based on the SER 
criterion an Eigen-beamformer is designed that exploits the 
knowledge of channel correlation available at the transmitter. 
The authors in 123 consider the same problem as in 1281 . 
yet quasi-orthogonal space time block coding (QSTBC) based 
beamforming is used instead of OSTBC. 

All works of ED-ISl consider the single-user multiple- 
input and multiple-output (MIMO) scenario. Recently, rank- 
two beamforming approaches have been independently pro¬ 
posed in l32l . l33l and l34l to enhance conventional single¬ 
group multicast beamforming in which multiple users are 
served on the same frequency resource subject to transmitted 
power and QoS constraints. A similar approach of rank-two 
beamforming is proposed for single-group multicasting using 
a relay network in 1351 . The concept of rank-two beamforming 
has later been extended in l3^ and 123 to solve the multi¬ 
group multicast beamforming problem in which beamforming 
is used to simultaneously deliver independent data services to 
users in predefined subscriber groups. 

By combining Alamouti coding with beamforming, rank- 
two beamforming approaches outperform the conventional 
rank-one approaches, however, the drawback associated with 
these rank-two beamforming approaches is that an optimal 
solution can only be obtained if the SDR solutiorQ exhibits 
a rank less than or equal to two, otherwise, an approximate 
solution is obtained in general. As discussed in 1331 . when the 

^With SDR solution we refer to the output obtained from interior point 
solvers such as CVX OH, ESI- 


rank of the SDR solution is greater than two, high dimensional 
(> 2) complex OSTBC can be applied instead of Alamouti 
coding to preserve the optimality of the beamforming solution, 
however at the expense of a reduced symbol rate associated 
with these OSTBC schemes BTl . 

The idea of combining beamforming with OSTBC in this 
work follows the general framework of 1221-123 in which 

rank-two beamformers are designed in combination with the 
application of Alamouti coding 12^ . In contrast to the rank- 
two beamformer designs in 123-123, we consider herein the 
downlink beamforming problem where each user is designed 
to be served by multiple beamformers combined with full-rate 
higher dimensional real-valued OSTBC. Real-valued OSTBC 
is employed because of its full (symbol)-rate property, thus the 
general rank approach proposed in the work achieves full-rate 
transmission as in the rank-one approaches of l24l . l25l and 
rank-two approaches of 123-123- 

In order to combine downlink beamforming with real-valued 
OSTBC, the effective channel vector of each user is adjusted to 
result in a real vector by applying a phase rotation procedure 
to which the optimal beamforming solution is proven to be 
invariant. Due to the orthogonality of the real-valued OSTBC, 
symbol by symbol detection can be performed at the receivers 
and the decoding complexity is not increased as compared 
to the conventional transmission that does not employ OS¬ 
TBC. The use of OSTBC results in multiple beamformers 
at each user and therefore multiplies the degrees of freedom 
in the beamformer design offering improved beamforming 
performance. Interestingly, the proposed beamformer design 
disembogues in the same SDR formulation as obtained in 
the conventional rank-one beamforming approaches of l24ll . 
l25ll and the rank-two beamforming approaches of 123-123, 
i.e., the beamforming problems after rank relaxation become 
identical. 

In the case that in the conventional rank-one downlink 
beamforming problem the rank of the SDR solution is greater 
than one, a rank reduction technique is applied to reduce 
the rank ll24ll . l25ll . l42l . Similarly, in the Alamouti coding 
based beamforming approaches of 123-123, rank reduction 
is applied if the SDR solution exhibits a rank greater than 
two. In our proposed real-valued OSTBC based beamforming 
approach, the SDR solution after the rank reduction procedure 
is proven to be optimal for the original problem if all ranks 
are no greater than eight. In the case that the SDR solution 
after rank reduction exhibits a rank greater than eight, random¬ 
ization techniques can be applied to compute an approximate 
solution Ca, 123- Eurthermore, we analytically prove that in 
our approach an optimal solution is always attainable if the 
number of additional shaping constraints does not exceed 79, 
whereas in the conventional rank-one approach in l24ll . l25ll 
and rank-two approach in 123-123, the maximal numbers 
of the shaping constraints are restricted to two and seven, 
respectively. Simulation results demonstrate the advantage of 
the proposed approach. 

The contribution of this paper can be summarized as fol¬ 
lows: 

• We address the problem of optimal QoS based downlink 
beamforming in the presence of a massive number of 
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arbitrary quadratic shaping constraints. 

• We increase the degrees of freedom in the beamformer 
design by combining optimal linear downlink beamform¬ 
ing with high dimensional real-valued OSTBC exploiting 
CSI knowledge at the transmitter. Our design can be con¬ 
sidered as a nontrivial full-rate extension of the Alamouti 
coding based rank-two beamforming framework of 1^ - 
Ezl to general rank beamforming supporting up to eight 
beamformers per user. 

• We analytically prove that in our approach an optimal 
beamforming solution can always be obtained when the 
number of additional shaping constraints does not exceed 
79. 

• Extensive simulation results demonstrate the effectiveness 
of the proposed downlink beamforming scheme in scenar¬ 
ios where the numbers of additional shaping constraints 
are extremely large. 

The remainder of the paper is organized as follows. Section 
II introduces the signal model and revises the conventional 
rank-one downlink beamforming problem. In Section III, the 
system model corresponding to the real-valued OSTBC based 
general rank beamforming approach is developed. Section 
IV formulates the optimal downlink beamforming problem 
involving real-valued OSTBC and provides the SDR solu¬ 
tion. Section V addresses the problem of computing optimal 
beamforming vectors from the SDR solution and provides a 
theoretic analysis regarding the optimality of the proposed 
downlink beamforming design. Simulation results are carried 
out in Section VI and conclusions are drawn in Section VII. 

Notation: E(-), || • || 2 , rank(-), Tr(-), diag{- ■ • }, 

and Z(-) denote statistical expectation, the Euclidean 
norm, the transpose, the Hermitian transpose, the rank of a 
matrix, the trace of a matrix, the diagonal matrix formed from 
the elements in the argument, the entry in the i-th row and 
j-th column of the matrix in the argument, and the argument 
of a complex number, respectively. Ik denotes the K x K 
identity matrix and X ^ 0 means that the matrix is a positive 
semidefinite matrix. denotes a sign in the set { >,<,=}■ 

II. Rank-one Beamforming 

We consider a cellular communication system where the 
serving base station equipped with an array of N antennas 
transmits independent information to M single-antenna re¬ 
ceivers. Let Si denote the information symbol for the t-th 
receiver with zero mean and unit variance. In conventional 
(rank-one) beamforming approaches of m-lU, ll24l . 1251 . the 
transmitter sends a superposition of signals {si}f£i for the 
different receivers using the respective iV x 1 beamforming 
vectors The received signal at the i-th single-antenna 

receiver is then given by a 

M 

yi= Siwfhi -f ^ (1) 

. T . ^1 

desired signal ^ ~ ^ 

interference plus noise 


variance af, respectively. The total transmitted power at the 

M 

base station equals ^ wf^w^. Based on ([T]), the SINK at the 
z-th receiver is derived as 


SINRi 4 


|wf h. 


M 

E 




( 2 ) 


Considering a QoS based beamforming design, we define ji as 
the minimum SINR requirement of the i-th user. Then the ex¬ 
tended downlink beamforming problem of minimizing the total 
transmitted power subject to minimum SINR constraints for 
each user and additional context specific shaping constraints 
can be formulated as ll24ll . ll2^ 


min 


s.t. 


M 

^ wf Wi 


M 


= 1,... ,M 

m—1 


(3a) 

(3b) 


M 

w^A/mWm >/ bi,yi=M+l,.. .,M+L{3c) 

m—1 

where Obl i represents a well-known reformulation of the QoS 
constraints with 



m = i, 
m ^ i, 


Vz, m = 1,..., M 
Vz, m = 1,..., M 


(4) 


and 


bi = ^i(Ti, >^ = >, Vz = l,...,M (5) 

and L additional quadratic shaping constraints are formulated 
in ( 1^ for appropriately chosen (as specified below) N x N 
Hermitian matrices Aj^jV/ = M + 1,..., M + L^ym = 
that are not necessarily positive definite, with 
corresponding thresholds 5j,V/ = M + -f L. Note 

that the shaping constraints in are not for information 
decoding purpose, while the first M constraints in Obl i are. 
Depending on the specific application under consideration, the 
additional shaping constraints in (O may take different forms 
(c.f. m, cs). Popular example applications which can be 
formulated under the framework of problem Q are described 
in subsections II.A and II.B. 


A. Positive Semidefinite Shaping Constraints 

In the context of cognitive radio networks. Aim — h;hp, 
where h/ denotes the channel vector between the base station 
and the Z-th primary user. In this case, with bi denoting the 
upper power threshold at the primary user and choosing = 
<, the Z-th general shaping constraint (l3^ takes the form0 

M 

y^ WmAim^ni < k- (6) 

m—1 


where h, and zz, are the A x 1 channel vector containing that if all matrices in S are positive semidefinite, then the 

the flat fading channel conditions and the receiver noise of problem j3) can be reformulated as a SOCP problem. 
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Thus the interference constraint (|6]l is used to guarantee that 
the power leakage to the primary users is below certain 
threshold a-m]. In the context of femtocell networks, h/ 
denotes the channel vector between the base station and the l- 
th concurrent user, and the shaping constraint (|6l) is designed to 
ensure that the power leakage to concurrent users in coexisting 
hierarchical networks is below certain threshold m-m. 

In the context of physical layer secrecy networks, in con¬ 
trast, h; denotes the channel vector between the base station 
and the Z-th eavesdropper, and the shaping constraint (|6l) is 
employed to enforce that the power leakage to eavesdroppers 
is below certain threshold ifTsll - lflTll . 

Similarly, in the context of energy harvesting networks, h/ 
denotes the channel vector between the base station and the 
Z-th charging terminal lfT8ll - ll20l . In this case 5/ denotes the 
minimum power threshold to be guaranteed at the charging 
terminal and = > is chosen. The shaping constraint (1^ 
can be rewritten as 

M 

m ^ hi (7) 

m—1 

with Aim — hihP for m = 1, M and Z = M +L 

to ensure efficient wireless charging. 

B. Indefinite Shaping Constraints 

Indefinite shaping constraints can be used to perform re¬ 
laxed nulling, as proposed in ||43]| . to reduce intercell inter¬ 
ference in multiuser downlink networks. Let h/ denote the 
channel vector from the base station of a given serving cell to 
a user of a different cell for which |he interference shall be 
limited. Defining Aim — hi = 0, and choosing 

/? as an appropriate interference threshold parameter lEI], the 
shaping constraint (1^ takes the form 

WmAim'W 

m ^ hi. (8) 

In this design the tolerable interference power induced by the 
Z-th user to the m-th user depends on the spatial signature h; 
of the co-channel user. Some other applications of indefinite 
shaping constraints can be found in BTl . Il22l . 


a, a. Since the SDR solution is not of rank one in general, 
rank reduction techniques are applied to obtain a solution to 
problem (|9]l with a reduced rank ll24l . ll25]l . see also Section 
V.A. However, in the case that a rank-one solution does not 
exist, an approximate solution can be computed from the SDR 
solution using, e.g., the popular randomization procedures of 
fM and ET). 


III. General Rank Beameorming 

The central idea of combining optimal downlink beam¬ 
forming with the concept of real-valued OSTBC proposed 
in this work follows the general framework of Ei-Ea 
in which rank-two beamformers are designed by combining 
beamforming with Alamouti coding and making use of CSI 
available at the transmitter. As compared to the rank-two 
approaches, we employ full-rate real-valued OSTBC to further 
increase the degrees of freedom in the beamformer design 
which grow linearly with the size of the code. Extending the 
rank-two beamforming approach to higher dimensional (>2) 
OSTBC has previously been considered as impractical due to 
the rate penalty associated with these codes ll^ . By applying 
real-valued OSTBC at the transmitter, multiple beamformers 
can be used to deliver the data stream to each user while 
maintaining the full-rate transmission property. 


A. Full-rate Real-valued OSTBC 

Let X (u) be a A X A real-valued OSTBC matrix given by 

ED 

K 

= (10) 

fc=i 

where A is the number of symbols per block, u = 
[ui,..., uk]"^ is an arbitrary A x 1 real vector and Cfc is a 
K X K real code coefficient matrix. Per definition the OSTBC 
matrix A(u) satisfies the orthogonality property 


C. Semidefinite Relaxation 

In this subsection we briefly revisit the SDR approach that is 
widely used to approximately solve the beamforming problem 
of form Q. The power minimization problem ([D is a quadrat- 
ically constrained quadratic programming (QCQP) problem 
which is NP-hard in general ll27l . Denote = Wiwf, 
problem Q can be rewritten as 


min 


s.t. 


M 


M 


Y, Tr(Aj^X,„) >ibi,Vl = l,...,M + L 

m—1 

X, ^ 0,rank(Xi) = 1, Vi = 1,..., M. 


(9) 


The SDR technique can be employed to solve the convex 
relaxation of problem (|9|l by removing the rank constraints 


A^(u)A(u) = X{u)X^{u) = ||u|| 2 Ik (H) 


which will be used in the following subsection. In this 
paper, we only consider real-valued OSTBC matrices with 
A= 1,2,4 or 8 which are the only possible sizes to achieve 
full rate ED . We note that, for a complex symbol vector u, the 
orthogonality property in (fTD can only be satisfied if A < 2 
Bril . Examples for real-valued OSTBC matrices are 


A([mi,U2,U3,M4]-' ) = 


T) A 

Ui 

U2 



_-U2 

Ui 

5 


Ui 

U2 

U3 

Ui 

-U2 

Ui 

—Ui 

U3 

-U3 

Ui 

Ui 

-U2 

— Ui 

-U3 

U2 

Ml 


( 12 ) 


, (13) 
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and X{[ui,... ,us]'^) = 


Ui 

U2 

Us 

Ui 

Us 

Ue 

Ur 

Us 

-U2 

Ui 

Ui 

-Us 

Ue 

-Us 

-Us 

Ur 

-U3 

— M 4 

Ui 

U2 

Ur 

Us 

-Us 

-Ue 

— U4 

Us 

-U2 

Ui 

Us 

-Ur 

Ue 

-Us 

-U5 

-Uq 

-Ur 

-Us 

Ui 

U2 

Us 

Ui 

-Uq 

U5 

-Us 

Ur 

-U2 

Ui 

— Ui 

Us 

-Ur 

Us 

U5 

-Uq 

-Us 

Ui 

Ui 

-U2 

-Us 

-Ur 

ue 

U5 

—Ui 

-Us 

U2 

Ui 


(14) 

B. System Model 

Denote = [s^i,..., SiK^ the K x \ complex symbol 
vector for the i-th user with K < N and K G {1,2,4, 8}, 
i.e., in correspondence with the dimension of the real-valued 
OSTBC matrices in (fT2li - (fT4li . In this work, we employ 
the real-valued OSTBC structure X{-) given in (flOl i on the 
complex symbol vector s^. Instead of weighting each symbol 
by a beamforming vector as in O, a code matrix T’(si) is 
transmitted for each user applying K beamformers of length 
N, denoted as Wa, ..., Wi/c . In this case, taking a slightly 
different perspective, each of the K beams can be regarded 
as a virtual antenna from which the OSTBC is transmitted. In 
our scenario we consider a block fading channel model where 
the channels remain constant over K time slots. The received 
signal Tjik at the i-th user in the fc-th time slot is given by 

M K 

Vik = + (15) 

m—1 k' — l 

where rnk is the noise of the i-th user in the fc-th time 
slot. In a compact matrix notation, the received signal vector 
yi=[2/a, • ■ •, at the i-th user within the transmission 

period of K time slots is given by 

M 

m—1 

M 

= T’(s,)Wfh,+ ^ T’(s^)W^h, + n, 

desiredsignal 

interference plus noise 

(16) 

where 

W, ^ iTe 11,2,4,8} (17) 

is the beamforming matrix, and ni=[nii, ■ • ■ ,'niK]'^- We as¬ 
sume that the noise vector at the i-th receiver is zero mean 
spatially and temporally white circular complex Gaussian with 
covariance matrix o-flR- The above system model can be 
reformulated in the following equivalent form M 

M 

y, = ^ T’(W^hj)s* -b n, 

m—1 

M 

= T’(Wfh,)s,-f ^ T'(W^h,)s„+n, (18) 

m = l,771^2 


where 


- A 

Yi = 

[yn, -ya, • ■ •, , 

(19) 


M 


T A 

u = 


(20) 


m—l,m^i 


- A 

= 

[^zl 7 '^z27 • ■ • 7 

(21) 


In order to implement full-rate transmission and symbol-wise 
decoding for each user, the code matrix X{yV^h.i) has to 
exhibit the orthogonality property (fTTl i. This however requires 
that the virtual channel vector becomes real-valued, 

i.e., the condition 

wfh, gM^, Vi = l,...,M (22) 

holds. We remark that in general is not real-valued for 

m i, and thus T’(W^hi) does not necessarily satisfy the 
orthogonal property in (fTTI) . Note that our proposed scheme 
is associated with a slight increase in the signaling overhead. 
This is due to the fact that the receivers have to know their 
individual composite channels for the decoding, i.e., the ith 
user requires the knowledge of Wf^h^. This signaling over¬ 
head is also required in the rank-one and rank-two schemes, 
however, the amount of signaling linearly grows with the code 
dimension of K. In this sense, the increase in the signaling 
overhead for our proposed scheme only applies to cases where 
there is a clear benefit in using general rank beamforming. 
In the following, we derive explicit expressions for the post 
detection SINR of the symbols received at the destinations 
under the assumption that condition (l22li is satisfied and that 
signal user detection is applied at the receivers. Towards this 
aim, we introduce the following useful lemma. 

Lemma 1. Assume that ip and w are a real and a complex 
vector both with the dimension K x 1, respectively. Let 
$ = X^{ip)X{(jj)X^{(jj)X{ip) where X{-) is a K X K real¬ 
valued OSTBC structure that fulfils SB- Then 

= Il'0ll2ll‘^ll2 Vfc = l,...,iT (23) 

where is the k-th diagonal element of the matrix 

Proof: Let X{uj) = Oi +jCt 2 where Oi and CI 2 are real 
orthogonal matrices from the definition of X{u>). Then 

X{uj)X^{uj) = (Hill}’ -f ) 

= \\u:\\llR+j{n2n'^ (24) 

Hence 

# = \\uj\\lx^{ip)x(ip)+jx^{ip){n2n{-nin^)x(ip) 

= (25) 

Since # is a Hermitian matrix and X{ip) is a real matrix, 
X^{ip){Cl 2 Cl^ — Cliri^)X(ip) is a skew symmetric matrix, 
i.e., its elements on the main diagonal are zero. Then the 
equation (|2^ holds. ■ 

For an orthogonal matrix XfW^h.i), i.e., with Wf^hj 
satisfying (l22T i. the transmitted symbol vector can be equalized 
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as 


l|Wfh,| 


A’"(Wfh,)y, 


— Si + 


1 


l|Wfh,||i 


A’"(Wfh,)(iz+n.). (26) 


Based on (l26l l. the covariance matrix of the received interfer¬ 
ence contained in is given by 

R?’ = ^ 


M 

E ^-(WfR)^(W-R 
X^{Wgh.)XiWfh,)] 


(27) 


and the covariance matrix of the noise in is given by 


r: 


(N) 


1 


l|Wfh,||4 


(Wf h,)E{n,nf }A’(Wf R) 


l|wfh,||: 


rlif ■ 


(28) 


Substituting according to (l22li the real-valued vector ■»/> = 
Wf^hi and complex vector u: = W^hi in (l27l i and applying 
Lemma 1, the interference power of the *-th user in the fc-th 
time slot can be expressed as 

= ||wfh.||j SllW.^‘‘-llf (29) 

With (|28]) and (|2^ . the post detection SINR corresponding to 
symbol Sik is given by 

E{sifcs*fe} 


SINR(s,fc) ^ 




|wfh,||l 


M 


(30) 


E ||W^R||i + af 

m = l,7717^2 


We note that the expression in ( |30] | is independent of the time 
index k and hence for the i-th user the post detection SINR is 
identical for all symbols in the OSTBC block. Thus we denote 
SINRi = SlNR(sife). For simplicity of presentation, the SINR 
constraints in the general rank approach can be written in a 
similar form as in the rank-one beamforming approach of Obl i. 
i.e., 


M 

E Tr(A,^W„W^) Vi = l,...,M (31) 

m—1 

where Aim is defined in (IHi. 

Since each symbol appears only once in each row of 
the code matrix T’(si), c.f. (fT2T i- (fT4l i. the transmitted power 
towards the i-th user in the fc-th time slot can be computed as 

P,k = E{efA(s,)WfW,A^(s,)e4 

= Tr(W,E{A^(s,)efeefA(s,)}Wf) 

= Tr(W,Wf) (32) 


power towards the i-th user in each time slot. Then the total 
transmitted power in each time slot amounts to 

M M 

E^2=ETr(W.Wf). (33) 

i=l i=l 

With multiple beamformers designed for each user instead of 
a single one, the additional shaping constraints in (1^ can be 
expressed as 

M K 

EE 

m—1 k—1 

M 

E Tr(W^Ai™W™) 

m=l 

M 

E Tr(Aj^W^W^) bi 

m—1 

VI = M + + L. (34) 

IV. The Power Minimization Problem 

In this section, we consider the problem of minimizing 
the total transmitted power per time slot subject to SINR 
constraints at each user and additional shaping constraints on 
the beamformers. Taking into account that according to (l22l) 
the virtual channel vector must be real-valued in order 

to satisfy the orthogonality property for simple decoding, the 
optimization problem is formulated in the following form 

M 

min VTr(W,Wf) (35a) 

M 

s.t. E Tr(A,^W„W") k 

m—1 

VI = 1,...,M + L (35b) 

Wf h, e , Vi = 1,..., Af. (35c) 

We remark that as a special case, the Alamouti code can 
be employed in our proposed scheme without the need of 
imposing the constraint (I35cl i. since the Alamouti code satisfies 
the orthogonality property of (fTTI) for an arbitrary complex 
vectors u while achieving full rate. In this case the proposed 
scheme becomes similar to the rank-two schemes proposed in 

El-El. 

A. Phase Rotation Invariance Property 

To solve the problem dTSl ). we first consider a relaxed 
problem of EH* by removing the constraints ( I35cl i. Let W* = 
[w”,...,w^] denote an optimal solution of the relaxed 
optimization problem of dTSl) in which constraint (I35cl) is 
omitted. Then we can perform the phase rotation on 
according to 

Wf 4 W*0, (36) 


where is the fc-th column of the N x N identity matrix. 
Similarly we observe that the transmitted power Pik is identi¬ 
cal in all K time slots. Let Pi = Pik represent the transmitted 


where the diagonal matrix &i is given by 
&i = diag{ exp(jZ(wf^hf),.. .,exp(jZ(w*f hf)}. (37) 
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Since satisfies all the constraints in ( l35T l, including 

constraint ( I35cb . it is a feasible solution to the unrelaxed 
problem ( |35] |. As the total transmitted power associated with 
is the same as that associated with the optimal 
solution we conclude that is an optimal 

solution to the original problem (l35T l. In other words, relaxing 
the real-valued requirements expressed in constraints (I35a 
in the beamforming problem ( |35] | results in an equivalent 
problem. An optimal solution for the original problem can 
always be computed from the solution of the relaxed problem 
by applying the phase rotation proposed in (l36T l. Therefore, 
we can, without loss of generality, omit constraint (I35cl l when 
solving (l35] l. 


B. SDR Approach 

Let us define the variable transformation 

X, 4 W,Wf. (38) 

Then substituting in (l35l l and adding the constraints 


X, ^ 0 (39a) 

rank(X,)<A:, Vi = l,...,M (39b) 

to ensure that the transformation (l38T l exists, problem (l35l l 
(with relaxed constraint OScb l converts to a rank constrained 
semidefinite program 

M 

^Tr(X,) (40a) 

M 

^Tr(Ai™XJ>i6i,VZ = l,... ,M+L (40b) 

m=l 

Xi ^ 0, (40c) 

rank(Xi)<A:, (40d) 

We remark that problem (l40b is identical to problem (|9]l 

except for the rank constraint. While in the latter problem 
the optimization variable X^ is restricted to the set of rank- 
one matrices, in our proposed formulation (l40b the rank of the 
matrix must not exceed K. This shows that the feasible set 
of our proposed beamforming approach is greater than that of 
the conventional one. 

Since the rank constraints in (l40b are non-convex, we 
employ the SDR approach a to obtain a relaxed convex 
optimization problem in which the rank constraints in (I40db 
are omitted, 

M 

min VTr(XA (41a) 

M 

s.t. ^Tr(Ai„X„)>i6i,V/ = l,...,M+L (41b) 

m—1 

X, ^0, = (41c) 


min 


s.t. 


problem of (1411 which has the following form [|24l 


max 

M + L 


s.t. 


where 


M+L 






M+L 


Zi =1- ^ piAii ^0 V* = 1, 


1=1 


pi>*0 Vl = l,...,M + L 






r> 

unrestricted 

< 




if is > 

if is = . 

if is < 


,M 

(42) 

(43) 


Note that, according to our previous observation, problem (HTb 
is identical to the SDR of the rank-one beamforming problem 
dUl. Therefore, the complexity at the transmitter side does not 
differ from that of the rank-one and rank-two schemes. This 
is due to the observation that the computational complexity of 
the proposed general rank approach mainly consists in solving 
(Eil, which is the same as in the rank-one and rank-two 
approaches. Problem (HTb belongs to the class of separable 
semidefinite programming (SDP) problems ll24l . ll25l and can 
be solved efficiently using solvers such as CVX EH, iol. 
Denote {X*})^]^ as an optimal solution to the problem (HTb . 
Then we can apply the rank reduction algorithm proposed in 
ll24l and ll25l with the input {X*}^]^ to reduce the rank of the 
optimal solution. A detailed description of the rank reduction 
procedure for general rank matrices is provided in Section V. 


V. Beamforming matrices generation 

In this section, we derive a sufficient condition on the max¬ 
imum number of shaping constraints under which a solution 
to (35) can always be obtained from the SDR solution. In this 
context, we adapt the rank reduction algorithm of ll24l . Il25l 
with modified stopping criteria for its application in general 
rank beamforming. Then we address the issue of determining 
the smallest code dimension K for all downlink users based 
on the output of the rank reduction procedure. In the case that 
a SDR solution after the rank reduction procedure has a rank 
greater than eight, a randomization procedure is proposed to 
obtain a suboptimal solution to the problem (l35l l. 


A. Rank Reduction Procedure 

The rank reduction procedure for general separable SDP 
of form @ has been proposed in ll24l . Il25l . By employing 
a modified stopping criteria, the rank reduction procedure is 
applied to compute a rank reduced solution from any optimal 
solution of (HTl i. 

Let {X*}^]^ denote a solution of the SDR problem dTIT i. 
The rank reduction algorithm successively reduces the rank 
of {X*}^j^ as follows. Introducing the factorization X* = 
QiQf where rank(X*) = rank(Qi) = Ki. Starting from the 


For later reference, we also provide the Lagrange dual the constraint is omitted. 
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given solution, the algorithm solves the following homoge¬ 
neous system of equations corresponding to ( |41bb 

M 

y] Tr(Q^Ai™Q„A^) = 0 I € {I,..., M + L} (44) 

m—1 

where G t^^K^xKm represents an unknown arbitrary 

Hermitian matrix. The number of real unknowns in (l44l i equals 

M 

^ rank^(X*), whereas the number of equations in (l44l i is 

2—1 _ 

M + L. Hence (144b must admit a nontrivial solution when the 
following inequality ll24l . Il25l 

M M 

Y, rank^ {X*) = Y k? < M + L (45) 

2=1 2 = 2 

is violated. A solution {X*}^^ that exhibits a reduced rank 
can then be computed as 

X* =Q,(I-^A0Qf (46) 

^max 

where (5max is the largest eigenvalue of all the matrices 
of {Ai}^r that satisfy (l44b . In ll24l . the above steps are 
repeated by assigning {X*}^^ to be a new input of the 
algorithm until the inequality (l43b is fulfilled. However, in 
our approach the rank reduction procedure is stopped after 
the maximum number of iterations max_iter is reached. The 
modified stopping criteria ensures that the ranks of {X*}^]^ 
cannot be further reduced when max_iter is set as follows 

M 

max_iter = ^rank(Xn-M, (47) 

i=l 

while the ranks still could be reducible if the stopping criteria 
in ll24l is employed. Note that for each iteration it requires 
to solve a homogeneous system of linear equations, matrix 
decomposition and singular value decomposition. But com¬ 
pared with the optimization problem (l4Tb . the operation cost 
of the iterations is comparatively small. The rank reduction 
procedure is summarized in Algorithm 1. 

Next we derive conditions on the number of additional 
shaping constraints and the code dimension K of the real¬ 
valued OSTBC for which optimal beamforming solution can 
always be obtained. These conditions are stated by the follow¬ 
ing lemma. 

Lemma 2. Assume that the relaxed problem (EB and its dual 
are solvabl^ and that the condition 

T < (A -f 1)2 - 2 (48) 

is satisfied, then there always exists an optimal solution X* 
for problem (ITO with rank(X*) < K for all i = 1,. .., M. 

Proof: We follow a similar line of argument as in ll24l and 
prove Lemma 2 by contradiction. Assume that (l48b is satisfied 

^“solvable” means that a bounded optimal value of the optimization 
problem can be obtained El. 


Algorithm 1 Rank reduction procedure 

Input {X*}^]^ an optimal solution to the problem (HTb . 

max_iter maximum number of iterations; 

Output {X*}^j^ such that the rank of any of the matrices 
{X*}ffi cannot be further reduced; 
while Number of iterations < max_iter do 
Decompose X* = QiQf Vi = 1,..., M\ 

Find a non-zero solution of the equation (l44b : 
if (l44b does not admit a nontrivial solution then 
break 
else 

Let bmax = where Su is 

l<i<M 

the /-th eigenvalue oF A^; 

Set X* = Q,(I - 3;^ A,)Qf Vz = l,...,M; 

end if 
end while 


and there exists a matrix X* with rank(X*) > K for some j 
such that the matrices 0 {X*}ff^ satisfy (fel l. Then 

(a) (b) 

Y rank2(X*) > M - 1 + {K + if > M + L (49) 

where strict equality holds in “(a)” if and only if there are 
M—1 rank-one matrices in {X*}ff^ and the last matrix has 
rank A-l-1, and the strict inequality in “(&)” follows from ( |48] |. 
The inequality (l49l l however contradicts our assumption that 
(|431 i is fulfilled. Hence all the matrices {X*}ffi must have 
ranks less than or equal to K. We conclude that the maximum 
number of additional shaping constraints L for which a rank 
less than or equal to K can be obtained is given by 

L= (K + lf -2. (50) 

■ 

Lemma 2 indicates that we can always find an optimal 
solution to problem dTSl l by using the SDR approach and the 
rank reduction procedure described in Algorithm 1 if condition 
(l50l l is satisfied. From (|50l), we can calculate the maximum 
numbers of additional shaping constraints for different choices 
of A G {1,2,4,8} as listed in Table |T] such that an optimal 
solution to problem (l35l l can always be obtained. We observe 
from Table U that our proposed scheme can accommodate a 
maximum number of 79 additional shaping constraints which 
corresponds to the choice of the code dimension A = 8. 

Since a smaller code size of the real-valued OSTBC matrix 
results in a shorter decoding latency at the receiver side, 
we seek to obtain the smallest value of A for all downlink 
users based on the output of the rank reduction procedure in 
Algorithm 1. If the updated {X*}ff^ after the rank reduction 
procedure satisfies rank(X*) < 8 for all i = 1,..., M, then 

®We observe that none of the matrices {X * }ff , are zero matrices, as 
otherwise at least one of the SINR constraints in 13 iF would be violated due to 
the positive semidefiniteness of and the definition of 

in 0- Hence all the matrices {X* must have a rank greater than or equal 
to one. 
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Number of beamformers 
per user K 

Number of additional 
shaping constraints 

1 

2 

2 

7 

4 

23 

8 

79 


TABLE I: Number of additional shaping constraints 


the smallest number K is chosen from K € {1,2,4, 8} such 
that 

A:>rank(X*) Vii = l,...,M. (51) 

The corresponding beamforming matrices are then obtained as 
= [Qi; OArx(i<'-rank(X*)] (52) 

where X* = Q,Qf with Q, e 

B. General Rank Randomization Procedure 

In the case that dril l is feasible, ( |48] | is violated and if at least 
one of the matrices in after rank reduction, exhibits a 

rank greater than eight, the following randomization technique 
which involves a power control problem could be applied 
to generate a feasible but generally suboptimal beamforming 
solution for problem dTSl l. Note that the randomization proce¬ 
dure may lead to an infeasible solution if the power control 
problem is infeasible. Note that the randomization procedure 
may not find a feasible suboptimal solution if the power 
control problem for each randomization sample is infeasible. 

Let us decompose the matrices {X*}^^ into 

X* = U,S,Uf. (53) 

The corresponding beamforming matrices are then 

randomly generated according to 

w, 4 [w,i, w, 2 ,..., w,8] = (54) 

where is the N x8 matrix whose elements are drawn from 
an i.i.d. complex circular Gaussian distribution with zero mean 
and unit variance. Note that the instances of the beamforming 
matrices generated in (|54] | are generally not feasible 

for problem dTSl l. In order to compute a feasible solution with 
spatial characteristics corresponding to {Wi}f£^, a power 
control problem involving linear programming is solved. The 
randomization procedure is summarized in Algorithm 2. 

With the rank reduction procedure in Algorithm 1 and the 
randomization procedure in Algorithm 2, a solution to problem 
(iTSl l can be computed following the procedure summarized in 
Algorithm 3. 

VI. Simulations 

Four simulation examples are provided to demonstrate the 
performance of our proposed downlink beamforming scheme 
with a large number of additional shaping constraints of 
different types. We assume that the base station is equipped 
with a uniform linear array (ULA), and the transmit antennas 
are spaced half wavelength apart. According to our system 
model, accurate CSI for all users and terminals is available 


Algorithm 2 Randomization procedure 

Input with ranks greater than 8 for some i, 

Aland number of iterations, 

Popt optimal value of the power control problem; 
Output The scaled beamforming matrices 
of the problem dril l: 

Set K = 8, Popt = +Inf; 
for fc = 1 to A{-and do 

Obtain according to ( l54l i: 

Solve the power control problem; 
if The optimal value is less than Popt then 
Set Popt to be equal to the optimal value and 
store the scaled beamforming matrices 
else 

Discard the matrices 

end if 
end for 


Algorithm 3 Summary 

Input {X*}^j^ an optimal solution to the problem dril l: 
Output {W*}^j beamforming matrices of the problem 
dS), K number of beamformers per user; 
if rank(X*) > 1 for some i then 

Apply Algorithm 1 to obtain the rank-reduced 
matrices {X*}^;^; 
end if 

if rank(X*) < 8 Vi = 1,..., M then 

Choose K to be the smallest number out of {1, 2,4, 8} 
such that rank(X* ) < K Vi = 1,..., M; 

Decompose {X*}^;^ to obtain using (15^ : 

else 

Apply Algorithm 2 to obtain suboptimal 
beamforming matrices 

end if 

Rotate matrices if necessary according to d22li . 


at the transmitter. The noise powers of the downlink users in 
all examples are assumed to be = 0.1 for i = 

We also declare that rank(Xj^) = ^ if the (^-|-l)-th largest 
eigenvalue is smaller than 0.01% of the sum of all eigenvalues. 

A. Example 1 

In the first example, we consider the design of downlink 
beamformers with external wireless charging terminals. Con¬ 
sidering a line-of-sight transmission scenario, we assume that 
three downlink users (M = 3) connected to the base station 
are located at directions 9i = —5°, 9^ = 10° and 9^ = 25° 
relative to the array broadside of the serving base station. 
The number of antennas at the base station is 12 {N = 12). 
We assume that there are 22 charging terminals, which are 
centered around 

I ? 4.....14 = [ - 80°, -75°, -70°, -65°, -60°, -55°, 

-45°,-35°,-25°,-8°,-2°] (55) 
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Fig. 1: The ranks of the matrices XJ (left bar), XJ (middle 
bar), Xg (right bar) after the rank reduction procedure. 

and 

t?i5,....25 = [12°, 18°, 35°, 45°, 50°, 55°, 

60°,65°,70°,75°,80°] (56) 

relative to the serving base station of the cell under considera¬ 
tion. For all downlink users and charging terminals, the spatial 
signatures are modeled as 

h(6l„) = 

me {1,...,25}, (57) 

i.e., the path loss of all downlink users and charging terminals 
is identical ll2^ . To make our simulation results more mean¬ 
ingful, we randomly vary the locations of the downlink users 
and the charging terminals in different Monte-Carlo runs, i.e., 
the angles of departure at the base station are simulated as 

dm =i}m + A9m m G {1,..., 25}, (58) 

where the random variations A9m drawn from a uniformly 
distributed within the interval [—0.25°, 0.25°]. We use the 
additional shaping constraints in (I41bl i to ensure predefined 
charging power levels at the l-th charging terminal in each 
time slot where Aim = h( 0 i)h^( 0 i) for m = 1 ,...,M and 
I = 4,..., 25. We set the minimum power threshold bi to be 
5dB for each charging terminal and = >. The SINK targets 
7 i at the individual downlink users are varied between OdB and 
lOdB. The simulation results are averaged over 300 Monte- 
Carlo runs. In each run, the number of randomization instances 
is set to 300 for all approaches if necessary. 

The ranks of the solution matrices of the relaxed problem 
(l4Tli after the rank reduction procedure are plotted in Fig. [T] 
According to the reduced rank property provided in Table |II| 
the code dimension is selected as K = 4. It can be analytically 
proven from a power scaling argument that problem (HTt , in 
the case of power charging constraints, is always feasible for 
all approaches. In Fig.|2]we display the total transmitted power 
per time slot at the base station versus the SINR for different 
approaches. As shown in Fig. |2| the proposed general rank 



Fig. 2; Transmitted power per time slot at the base station. 

beamforming approach outperforms the competing approaches 
in terms of transmitted power. In the low SINR region, the 
gap between the rank-one and rank-two approaches and the 
proposed approach is large, because as shown in Fig. [T] most 
of the solution matrices are of high rank (> 2 ) and thus the 
suboptimal randomization approximation is performed in the 
rank-one and rank-two approaches. In the high SINR region, 
the gaps between different approaches decrease because as 
shown in Fig. [T] the percentage of rank-one solution matrices 
increases which results in an increased number of optimal 
solutions for all approaches. 

B. Example 2 

In the second example, we consider the downlink beam- 
former design according to problem dril l for beam pattern (BP) 
with smooth and flat sidelobes to reduce the interference to 
co-channel users. We assume that in our simulation scenario 
the base station consists of 18 antennas (N = 18). In this 
simulation the locations of three downlink users (M = 3) are 
the same as in the previous example, i.e., 9\ = —5°, 02 = 10° 
and 03 = 25°. The SINR thresholds of downlink users are set 
to = lOdB. Moreover, we assume that nineteen co-channel 
users connected to a neighboring base station are present in 
the scenario, which are located at 

Mi....,i 9 = [-89.375°, -80°, -70.625°, -61.25°, -51.875°, 
- 42.5°, -33.125°, -23.75°, -14.375°, 2°, 3°, 17°, 18°, 
34.375°, 43.75°, 53.125°, 62.5°, 71.875°, 81.25°]. (59) 

The channel propagation model is the same as defined in dSTl) . 
The interference power at the direction fii relative to the base 
station in each time slot can be written as 

M 

Tr(A(j+ 3 )„X,„) (60) 

m—1 

where = h.{pj)h.^{pj) for m = 1 ,...,M and j = 

1,..., 19. In our beamformer design, the interference power 
is upper bounded by 6^+3 = 0.1 and ^_,+3 = < for j = 
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X) 

X) 

x.^ 

Original rank in <41i 

14 

15 

15 

Reduced rank 

2 

3 

4 


TABLE II: Rank property before and after applying rank 
reduction algorithm. 


1,19. In addition to these constraints, we guarantee that the 
interference power at the direction attains a local minimum 
value by adding interference derivative constraints, i.e., the 
interference in the vicinity of the constraint directions remains 
approximately constant if 




and 




d^j 

J 19} 

where the threshold is set to Ca = 10 “®, 


d-K 


> 0 


( 61 ) 


dfij 

and 

dp? 


M 


— j G {20,..., 38} (62) 


M 


= ^ Tr(A(,+ 3 )^X„), j e {58,..., 76} (63) 


are satisfied, for m = 1,..., M, 

+ 3 = 20,38 

j=39,...,57 


Ki+3)» 


, = < 


dfi 

dfij 


‘3^ d^j 




'^3'^ d^j 

3) 


dfij dfi:- ’ 


('i+3 - 


, j=20,...,38 

Ca, } = 39, ...,57 
0, j=58,...,76 

<, j=20,...,38 
>j+ 3 =<(>, j=39,...,57 
>, j=58,...,76 


j = 58,...,76 

(64) 

(65) 


( 66 ) 


with j = j mod 19, i.e., the remainder of j divided by 19. 
The received sum power at direction 0 relative to the base 
station, referred to as the sum BP, is defined as 

M 

Y, l|h(d)w:,||^ (67) 

m—l 

where is the rank reduced solution given in (l52l i. 

The BPs are presented in Fig. |2 and the rank properties of 
the solution matrices are provided in Table |II] Note that there 
is a total number of 76 additional shaping constraints in this 
simulation. According to Lemma 2, we can find an optimal 
solution to the optimization problem (HTt with the rank less 
than or equal to 8 by using the rank reduction procedure. Based 
on the results in Table HU we select the code dimension K = A. 

As shown in Fig. |3l the proposed approach is capable of 
coping with a large number of additional shaping constraints. 



Fig. 3: User BPs and sum BP with smoothed and suppressed 
sidelobes. 


Furthermore, as listed in Table [III the ranks of the solution 
matrices have been significantly reduced which demonstrate 
the effectiveness of the rank reduction procedure. 


C. Example 3 

The same scenario as in Example 2 is considered to 
perform a comparison between our proposed approach with 
the conventional rank-one and rank-two approaches. All 
location parameters remain unchanged. Furthermore, we 
assume that all angles of departures are also subject to 
variations in different Monte-Carlo runs, which are defined 
in the same way as in Example 1. The required SINRs 
at the downlink users are uniformly varied between OdB 
and 5dB. The results are averaged over 300 independent 
Monte-Carlo runs and the number of randomization instances 
in each run is set to 100 for all approaches if necessary. The 
feasibility percentage of all approaches is displayed in Fig. 
ID From Fig. |4] we observe that the proposed approach is 
always feasible for different SINR thresholds. In contrast to 
this, the feasibility of the rank-one and rank-two approaches 
decreases with increasing SINR thresholds. This demonstrates 
that our proposed approach has a wider feasibility range 
compared to existing approaches. The ranks of the solution 
matrices of the relaxed problem drill after the rank reduction 
procedure are plotted in Fig. |5] As shown in Fig. 4, when 
7 i < 3dB, all three approaches are feasible. This is due 
to the fact that in this case, as shown in Fig. 5, rank-one 
solutions are obtained for all approaches. In other words, 
optimal solutions are obtained for all approaches and thus 
the performance obtained from all approaches is identical. 
Therefore, when 7 ^ < 3dB, the code dimension for our 
proposed method is chosen as AT = 1. In contrast to this 
when 7 i > 3dB we observe from Fig. 5 that the rank of the 
optimal solutions takes different values in the range between 
one and five. Thus in contrast to the rank-one and rank-two 
beamforming approaches if a rank larger than two is obtained, 
our proposed approach retains the optimality property and 
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0(degree) 


Fig. 4: The feasibility percentage of all approaches. 


Fig. 6 : User BPs and sum BP of downlink beamforming 
problem subject to interference power constraints. 



0 1 2 3 4 5 


SINR(dB) 

Fig. 5: The ranks of the matrices (left bar), XJ (middle 
bar), Xg (right bar) after the rank reduction procedure. 


yields feasible solutions while the competing approaches yield 
suboptimal solutions or even become infeasible for 7 ^ > 3dB. 


the interference power to the coexisting users by the following 
constraints 

Tr(h(0,)h(0,)^X.) < /3||h(0,)||iTr(XO 
Vi = 1,2,3, Vj = 4, ...,25 ( 68 ) 

where /? <C 1 is an interference constraint parameter. The 
above constraints can be reformulated into the form of (I41bl i 
where, for m, m = 1,..., 3, 

* [/3||h(6'j)||il-h(6»j)h(6»j)^, m = m 

m ^ m 

(69) 

(70) 

^n +3 =>, n = 1, ..., 66 ; j = 4, ...,25. (71) 

We note that the matrix Aim is either zero or indefinite for 
I = 4, ...,69,m = 1,...,3 and there is a total number of 
66 additional shaping constraints in this simulation. In the 
simulation, (3 is chosen to be 0.5%. As shown in Fig. | 6 l the 
interference power at the locations of the coexisting users is 
limited to a reasonable level. 


^(22(m-l)+j)m = 


0 , 


bn+3 — 0 , 


D. Example 4 

The aim of the fourth example is to demonstrate the 
interference power suppression at each co-channel user to a 
fraction of its maximum value. In this example the concept of 
relaxed nulling is used to formulate the additional (indefinite) 
shaping constraints for interference power limitation Il2n . The 
base station under consideration is equipped with a ULA of 
15 antennas that are spaced half wavelength apart {N = 15). 
Three downlink users served by the base station are located 
at 9i = —15°, 02 = 5° and 0^ = 25° relative to the base 
station. We assume that twenty two co-channel users served 
by neighboring base stations are present in our scenario which 
are located at the same position as in dSST l and (|56] |. We set the 
SINK thresholds to the same value as in Example 2. Similarly, 
the spatial signatures are modeled according to (1571) . We limit 


VII. Conclusion and outlook 

In the paper, we propose a general rank beamforming 
approach for the multiuser downlink beamforming problem 
with additional shaping constraints. The general rank approach 
increases the degrees of freedom in the beamformer design by 
using high dimensional full-rate real-valued OSTBC. In our 
proposed approach, an optimal solution can be obtained when 
the ranks of all SDR solution matrices are less than or equal 
to eight after the rank reduction procedure. Moreover, in our 
scheme an optimal solution for the original problem can be 
found when the number of additional shaping constraints is 
less than or equal to 79. The range of applications for our 
proposed beamforming scheme is hence much wider than that 
of the conventional rank-one and rank-two approaches. Our 
proposed general rank beamforming framework exhibits an 
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underlying optimization problem structure that is similar to 
that of the conventional rank-one and rank-two beamforming 
approaches. This allows a simple extension of the approach to 
existing robust beamforming designs, e.g., in the practically 
important case of inaccurate CSI, maintaining the benefits of 
increased degrees of freedom available in the proposed general 
rank beamforming approach. All the results presented in this 
paper can be extended to QCQP problems with double-sided 
constraints considered in ll45l as an interesting topic of future 
research. 
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